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Abstract 

In this paper we present calculations of thermodynamic functions within Zhang's SO(5) quantum 
rotor theory of high-T c superconductivity. Using the spherical approach for the three-dimensional 
quantum rotors we derieved explicit analytical formulas for entropy and specific heat related to the 
lattice version of the SO (5) nonlinear quantum-sigma model. We present the temperature depen- 
dence of these quantities for various settings of relevant control parameters (quantum fluctuations, 
chemical potential). We find our results in overal qualitative agreement with basic thermodynamics 
of high-T c cuprates. 
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I. INTRODUCTION 



The theory unifying antiferromagnetism (AF) and superconductivity (SC) to describe 
global phase diagram of high-T c superconductors was recently proposed by Zhang.0 In this 
approach, based on symmetry principles, a three-dimensional order parameter (the stag- 
gered magnetization) describing the AF phase and a complex order parameter (with two 
real components), describing a spin singlet d-wave SC phase are grouped in five-component 
vector called a "superspin". The SO(3) symmetry of spin rotations (which is spontaneously 
broken in the AF phase) and the electro- magnetic SO (2) invariance (whose breaking defines 
SC phase) along with well-defined AF to SC and vice versa rotation operators form SO (5) 
symmetry. In the Zhang's theory both ordered phase arise once SO (5) is spontaneously 
broken and the competition between antiferromagnetism and superconductivity is related 
to direction of the "superspin" in the five-dimensional space. The low energy dynamics of 
the system is determined in terms of the Goldstone bosons and their interactions specified 
by the SO (5) symmetry. The kinetic energy of the system is that of a SO (5) rigid rotor 
and the system is described by a SO (5) non-linear quantum a model (NLQctM). The SO (5) 
quantum rotor model offers a Landau-Ginzburg-like (LG) approach for the high-T c prob- 
lem. However, it goes much beyond the traditional LG theory, since it captures dynamics. 
While the SO (5) symmetry was originally proposed in the context of an effective field-theory 
description of the high-T c superconductors, its prediction can also be tested within micro- 
scopic models.i'i'B&ifl For example, numerical evidence for approximate SO (5) symmetry 
of the Hubbard model came out from exact diagonalization of small sized clusters.! The 
global features of the phase diagram deduced from SO (5) theory based on a spherical quan- 
tum rotors! agree qualitatively with the general topology of the observed phase diagram 
of high-T c superconductors. The quantitative investigation of the quantum critical point 
scenario within the concept of the SO (5) group, e.g. the scaling of the contribution to the 
electrical resistivity due spin fluctuations showed linear resistivity dependence on tempera- 
ture for increasing quantum fluctuation - being a hallmark example of anomalous properties 
in cuprate materials.0 The systematic studies of magnetic properties of the SO (5) theory 
showed that the theory yields a qualitative scenario for the evolution of magnetic behavior, 
which is consistent with experiments.0 It qualitatively explains the results of experimen- 
tal measurements (notably the NMR relaxation rates) with correct predictions of behavior 
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of uniform spin susceptibility in high temperatures. Also the energy dependence of the 
momentum-integrated dynamic spin susceptibility show features, which are in qualitative 
agreement with experimental findings. 

Thermal fluctuations are pronounced in the high-T c superconductors due to number of 
reasons. The carrier density is rather small, the anisotropy is large and the critical temper- 
ature is high. It turns out that deviations from the mean-field behavior are present in the 
specific heat C at the superconducting transition temperature T c . In the mean-field BCS 
theory, a second order transition with a jump in a specific heat at T c takes place. In contrary, 
in most high-T c superconductors thermal fluctuations seem to restore common behavior. 

The aim of the present paper is to study quantitatively basic thermodynamic functions 
resulting from the SO (5) theory, thereby substantiating this theoretical framework. Our 
study may also provide a useful diagnostic tool for testing the basic principles of SO (5) 
theory by comparing the quantitative predictions {e.g., specific heat) with the outcome of 
the relevant experiments. 

The outline of the reminder of the paper is as follows. In Section II we begin by setting 
up the quantum SO (5) Hamiltonian and the corresponding Lagrangian. In Section III we 
find closed forms of various thermodynamic functions. We calculate free energy, entropy 
and specific heat. Finally, in Section IV we summarize the conclusions to be drawn from 
out work. 



II. HAMILTONIAN AND THE EFFECTIVE LAGRANGIAN 

We consider the low-energy Hamiltonian of superspins placed in the nodes of a discrete 
three dimensional simple cubic (3DSC) lattice, 

ZU i \1<V i<j 

-V(n l )-2^L] 5 - (1) 

i 

Indices % and j number lattice sites running from 1 to N - the total number of sites, while 
fi,u = 1,...,5 denote superspin rij = (ni, n 2 , n 3 , n 4 , n 5 ) i components (n^F,* = (n 2 ,n 3 ,n 4 ) i 
refers to antiferromagnetic and nsc,i = (^i,^s)j superconducting order, respectively). The 
superspin components are mutually commuting (according to Zhang's formulation) and their 
values are restricted by the rigidity constraint nf — 1. 
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The first part of the equation ([]]) is the kinetic energy of the system (being simply that 
of a SO (5) rigid rotor), where 

= n niPui ~ n uiPiii (2) 

are generators of Lie SO (5) algebra (expressed by total charge L} 5 , spin and so-called 
"7r" operators), p^i are momenta conjugated to respective superspin components: 

d 



Via 



[n^Pu] = iS^ (3) 

and parameter u measures the kinetic energy of the rotors (an analog of moment of inertia). 

The second part of the Hamiltonian is the inter-superspin interaction energy with J being 
the stiffness in the charge and spin channel. In the 3DSC lattice, J is nonvanishing for the 
nearest neighbors and its Fourier transform 

is simply 

J q / J = cos q x + cos q y + cos q z . (5) 
For convenience, we will further introduce the density of states 

p(0 = ^EHZ-JJJ), (6) 

q 



which for 3DS0i=l lattice reads 



1 /-™(i.2-0 dy 

TT^ Jmax(-l,-2-£) y/1 — y z 



xK 



\ 



0(3- kl), (7) 



where K (x) is the elliptic integral of the first kind and O (x) is the step function^. 

The last two parts of the equation ([!]) provide symmetry SO (5) breaking terms. In the 
result of their interplay, the system favours either the "easy plane" in the superconducting 
(SC) space (ni, n§), or an "easy sphere" in the antiferromagnetic (AF) space (712, n.3, 114). The 
first of two terms is defined as: 



v ( n = f E K + < + O > ( 8 ) 



with the anisotropy constant w, which positive value favours AF state. The second term 
contains a charge operator Lj 5 , whose expectation value yields the doping concentration and 
the chemical potential // (measured from the half- filling) , whose positive values favour SC 
state. 

The partition function Z = Tre~ H ' kaT is expressed using the functional integral in the 
Matsubara "imaginary time"r formulation! (0 < r < 1/ksT = j3, with T being the temper- 
ature). We obtain: 



n p*] / n 



Dpi 



2tt 



6(1- n 2 { ) 6 (m ■ Pi) x 



= Jl[[Dn t }6(l-^)e-iy^\ 

i 

with £ being the Lagrangian: 



(9) 



C(n) 



fdn sc \ . fdn AF \ 



2 2 , • / 

-4-u/i n sc + Aiujj, \ -^—n^ 



dr 



dn 5 
dr 



■m 



+ 



Jij^i ■ % - — ( n li + n li + n li) 
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The problem can be solved exactly in terms of the spherical mo defl To accommodate 
this we notice that the superspin rigidity constraint (nf = 1 ) implies that a weaker condition 
also holds, namely: 



JV 



(ii) 



i=l 



Therefore, the superspin components n, (r) must be treated as c-number fields, which 
satisfy the quantum periodic boundary condition rij (j3) = rij (0) and are taken as continuous 
variables, i.e., — oo < n.; (r) < oo, but constrained (on average, due to Eq. ([□])) to have 
unit length. 

This introduces the Lagrange multiplier A (r) adding an additional quadratic term (in rij 
fields) to the Lagrangian fllPl). Fourier transform n(k, ug) of the superspin components 



n,- T 



-| oo 

— £ £ n(k,^)e-^-- k -) 

P 1 ^ k £=-oo 



(12) 



introduces the Matsubara (Bose) frequencies u>e = 2tt£/P (£ = 0, ±1, ±2, 
Using the Eq. (|9]), the partition function can be written in the form: 



-N<t>{\) 



2ni 



(13) 



where the function <p (A) is defined as: 



0(A) 



^drA (r) - i In / n [^1 x 



x exp 



(14) 



The exact value of the partition function can be found in the thermodynamic limit (N — > 
oo), when the method of steepest descents is exact and the saddle point A (r) = Ao satisfies 
the condition: 

54>(X) 



8X(t) 



0. 



(15) 



A=A 



At criticality, corresponding order parameters susceptibilities becomes infinite and corre- 
sponding Lagrange multipliers are: 



A 



AF 



A 



sc 



1 w 
2 Jk=0 + 2 ' 

1 , 

- J k =o + 2xfJ> , 



(16) 



for AF and SC critical lines, respectively. Furthermore, using the spherical condition ( [TT]) 
and the values (p!6|), we finally arrive at the expression for the critical lines separating AF, 
SC and QD (quantum disordered) states (for more specific description these calculations 
see, Ref. ggj). 



III. THERMODYNAMIC FUNCTIONS 



A. Free energy 



The free energy is defined as / = - (/5A)" 1 InZ = (/?P <f) (A ). Using the formula (0) 
we obtain: 



/ 



-A + 



2pN 



[2A- J k + 



UUJo 



W 



+ 



(17) 



+ ^ E ln [ 2A - J k + u (u t + 2tfi) 2 } + -i_ ^ In [2A - J k + n (u t - 2tfif 
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After performing the summation over Matsubara's frequencies, we obtain the free energy: 



+ In 2 sinh 
+ In 2 sinh 



2A- 




u 


'2A- 




it 



2/i 



+ 2/i 



B. Entropy 

The entropy is defined as S = k B [3 2 ^. Using the formula (|1^) we obtain: 



S(P) 



k B 
2 

+ 
+ 



f 00 PiOdds 

J — oo ( 



(3 A (0 coth (0 - 2 In 2 sinh ^A (f) 



+ 



l3B_ (0 coth (0 - 2 In 2 sinh (0 
(0 coth (0 - 2 In 2 sinh ^B + (f) 



where: 



(18) 



(19) 



A® 
B- (0 

5+ (0 



'2A- Jf - w 



'2A 




u 


>2\ 





-2/i, 
+ 2/x. 



(20) 



The dependence of the entropy on temperature and chemical-potential is shown on Fig. [I]. 
Starting from T = 0, the entropy increases in any ordered phase (AF or SC) until reaching 
T c (or T/v). The further increase is slower, but saturation in higher temperatures is not 
observed. The absolute value of the entropy is lower for higher quantum fluctuation (see, 
Fig. |^). We find obtained results in qualitative agreement with experimentally measured 
properties of high-T c superconductors (e.g. for Bi2212 compound, see Ref. O). 
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k B T/J 



FIG. 1: Plot of the entropy S vs. chemical potential [i/J and temperature k B T/J for fixed uJ = 3 
and w/J=l. Solid lines indicate the projection of the [i — T phase diagram. 
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FIG. 2: Plot of the entropy S vs. temperature k B T/J for w/J = 1, fi/J = 0.2 and different values 
of uJ, as indicated in the inset. 

C. Specific heat 



The specific heat at constant volume is defined: 

d 2 



C = -k B (3 z 
dX 



— {pf) = -k B p 

'd 2 fdX +2 _d 2 f 



' df d 2 f 



dX 2 df3 d\d(3_ 

The derivative dX/d(5 can be found from the saddle-point condition (|15|) : 

df [A 09) , 0\ 



dX 



0. 



(21) 



(22) 



A=A 



s 



Explicitly, we obtain: 



The specific heat: 



dX 

d(3 



a 2 f 

dXd/3 
d 2 f 
d\ 2 



c = -k B p 2 
Using the formula (|18|) we obtain: 



df B 2 f d\ d 2 f ' 
2— 4- B — 4- B — 

d(3 p d(3 2 p d(3d\d(3 



(23) 



(24) 



x |3A 2 (0 sinh" 2 (0 + Bl (0 sinh" 2 (£) + B\ (£) sinh" 2 ^B + (£) J + 
P (0 ^ x 



4m cf/3 j-oo 

x | 3 sinh" 2 I A (0 + ^ sinh" 2 (0 + ^ sinh" 2 | S+ (0 } , (25) 
where A (£) , 5_ (£ ) and _B + (£ ) are defined by formulas (pO"D . 



C(Z) = \r-±, (26) 



and 



g-tf) ■ ,,-2/3 R m , g+(0 ■ h -2i3 R m \, 
_ smh _B_ K) + _ 8mh -B + ({)|/ 

The temperature dependence of the specific is presented in Fig. |3|. The low temperature 
behavior of C(T) may be approximated by C(T) ~ T 3 for /j,/J = and C(T) ~ T 2 ' 5 for 
/i/J = 0.44. For higher temperatures (but still below transition temperature) the linear 
behavior of the specific heat is observed. Reaching the critical temperature (T c or TV), the 
specific heat experiences a finite jump (implying the value a = for the critical exponent 
of the specific heat). For higher temperatures, saturation is observed. 
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FIG. 3: Specific heat vs. temperature fcgT/J for w/J = 1, uJ = 3 and various values of chemical- 
potential, as indicated on the figure. 

IV. SUMMARY AND FINAL REMARKS 

In conclusion, we have calculated entropy and specific heat dependence on temperature 
and various other parameters using the unified theory of antiferromagnetism and supercon- 
ductivity proposed for the high-T c cuprates by Zhang and based on the SO (5) symmetry 
between antiferromagnetic and superconducting states. The theory of yields a qualitative 
scenario for the evolution of thermodynamic functions behavior, which is consistent with 
experiments. Most experimental work on the specific heat in the high-T c superconductors 
have concentrated on yttrium compound (Y-123). Optimally doped Y-123 does not 
show a jump in the specific heat, but A-peak at T c . However, the shape for overdoped Y- 
123 is intermediate between a BCS step and a A-type transition. Furthermore, optimally 
doped Bi-2212 shows a symmetric anomaly (intermediate between A-peak and finite jump). 
Experimentally, the specific heat is not very sensitive to the critical exponent a and one 
can ascertain that |a| <C 1 for Y-123 compounds. However, the result of the present work 
(a = 0) agrees with the critical behavior of the 3D-XY model. Finally, checking the validity 
of basic principles of the SO (5) theory, by comparing parameters discussed here with rele- 
vant, obtained from calculations on microscopic models of high-T c superconductors is still 
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called for. 
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